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SUMMARY 

This  report  describes  a new  method  for  obtaining  the  low 
incidence,  linear  aerodynamic  data  for  a missile  from  very 
simple  free  flight  trials  on  a gas  gun  range.  The  only 
instrumentation  in  the  free  flight  model  is  two  flashing 
lights,  one  in  the  nose  and  one  in  the  tail.  The  images  of 
these  flashing  li^ts  are  recorded  by  three  ballistic 
cameras  and  the  positions  of  the  lights  can  then  be  cal- 
culated very  accurately  by  a least  squares  solution  of  the 
resulting  triangulation  problem.  The  velocity  of  the  centre 
of  gravity  and  the  attitude  of  the  vehicle  can  be  estimated 
from  the  position  measurements  for  the  nose  and  the  tail. 

A parameter  estimation  technique  using  the  output  error 
criterion  is  applied  to  the  attitude  and  velocity  data  to 
obtain  aerodynamic  terms.  The  method  has  been  used 
successfully  to  estimate  axial  force,  normal  force,  static 
stability  and  pitch  da^)ing  moment  for  the  test  vehicle. 
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1.  INTRODUCTION 

The  concept  of  a double  flashing  light  vehicle  grew  from  the  single  flashing 
light  tests  reported  in  references  1 and  2.  The  single  flashing  light  technique 
has  been  used  to  obtain  very  accurate  trajectories.  The  position  data  obtained 
with  the  single  flashing  light  were  so  accurate  that  it  seemed  likely  that  useful 
attitude  data  could  be  obtained  if  longitudinally  separated  lights  were  carried. 
The  experimental  method  is  described  in  detail  in  Section  2.  The  principal 
ingredients  are  two  flashing  lights,  which  are  carried  in  the  nose  and  tail  of 
the  missile  and  three  ground  based  ballistic  cameras  which  are  used  to  record 
the  flashes  along  the  vehicle  trajectory.  A set  of  four  reference  lights  is 
used  to  calibrate  the  camera  plates  and  the  position  of  each  flash  is  derived. 

The  nose  and  tail  positions  obtained  in  this  way  are  converted  to  centre  of 
gravity  trajectory  data  and  attitude  da^’a,  for  the  vehicle. 

The  first  flight  of  a vehicle  carrying  two  flashing  lights  was  reported  in 
reference  1.  The  vehicle  was  launched  at  zero  incidence  and  the  derived  values 
of  incidence  were  regarded  as  indicative  of  the  r.m.s.  noise  levels  and  bias 
errors  which  might  be  present  in  the  data  from  such  trials.  On  that  basis  this 
preliminary  trial  showed  that  attitude  angles  could  be  measured  with  r.m.s.  noise 
levels  of  less  than  one  degree  and  bias  errors  of  about  half  a degree.  However, 
it  was  expected  that  the  bias  errors  could  be  substantially  reduced  by  a more 
accurate  survey  of  the  range. 

The  unexpectedly  high  quality  of  the  trajectory  and  attitude  data  obtained 
from  the  first  trial  prompted  some  rethinking  on  the  uses  of  the  technique.  As 
consequence  the  second  vehicle  was  launched  at  incidence  with  the  aim  of 
:racting  as  much  aerodynamic  data  as  possible  from  the  trial.  The  data 
analysis  has  been  carried  out  using  a parameter  estimation  technique  which  relies 
on  an  output  error  criterion.  The  method  attempts  to  find  values  for  the 
parameters  characterising  the  mathematical  model  such  that  the  sum  of  the  squares 
of  the  differences  between  the  model  outputs  and  the  measured  values  for  attitude 
angles  and  centre  of  gravity  velocity  components  is  minimised.  The  implement- 
ation of  the  parameter  estimation  technique  is  outlined  in  Section  3. 

Unfortunately,  no  measurements  of  roll  orientation  were  made  on  this  trial, 
so  that  only  roll  independent  aerodynamic  terms  could  be  included  in  the  model. 
Thus  the  model  cannot  represent  non-linear  aerodynamic  terms  effectively. 

Although  minimal  non-linear  effects  are  included  in  the  mathematical  model  of 
Section  3,  they  are  not  well  determined  by  the  data  analysis.  Over  90  per  cent 
of  the  data  points  have  incidence  values  less  than  six  degrees.  Wind  tunnel 
data  from  reference  3 indicates  that  non-linear  effects  are  negligible  in  this 
incidence  range.  The  values  obtained  for  the  parameters  are  presented  in 
Section  4 and  comparisons  are  made  with  wind  tunnel  results.  The  accuracies 
obtained  are  surprising  in  view  of  the  quality  of  the  data;  r.m.s.  noise  levels 
vary  from  1%  in  the  static  aerodynamic  pitching  moment  derivative  iq>  to  7%  in 
the  dynamic  pitch  damping  moment  derivative..  The  results  obtained  are 
particularly  encouraging  because  the  r.m.s.  noise  level  in  the  measurements  of 
attitude  angles  is  between  three  and  four  times  the  noise  level  in  the  first 
experiment.  Hence  we  expect  that  results  from  future  trials  will  be  much 
iiqjroved. 

Section  5 presents  some  conclusions  which  are  drawn  from  the  general  success 
of  the  trials  technique  and  the  data  analysis.  The  principal  conclusion 
concerns  the  next  step  in  the  development  of  the  method.  The  method  is  restrict- 
ed, by  the  assumption  of  linear  aerodynamics,  to  treating  situations  where 
incidence  amplitudes  are  small.  In  order  to  extend  the  method  to  treat  non- 
linear aerodynamic  terms  adequately,  some  measure  of  roll  orientation  is  required, 
since  the  non-linear  contributions  to  the  aerodynamics  of  finned  vehicles  always 
depend  on  roll  angle.  Plans  are  presently  under  way  for  testing  methods  of 
measuring  roll  orientation. 
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2.  TRIAL 

The  test  vehicle  is  depicted  in  figure  1 and  physical  data  on  the  missile  are 
given  in  '^.ible  1.  It  was  launched  from  a 385  naa  bore  compressed  air  gun  at  a 

velocity  of  120  m/s  and  an  elevation  of  20°  on  a small  gas  gun  range.  The  lay- 
out of  the  range  is  shown  in  figure  2. 

The  launching  technique  attempts  to  induce  the  missile  incidence  to  oscillate. 
It  is  arranged  in  the  sabot  at  an  angle  of  about  eight  degrees,  and  as  it  emerges 
from  the  sabot  air  flows  only  over  the  forebody.  The  general  effect  is 
destabilising  and  the  incidence  increases  until  flow  becomes  established  over  the 
tail.  In  this  way  initial  incidences  up  to  20°  can  be  obtained.  This  provides 
an  adequate  pitching  motion  for  later  analysis.  The  nose  and  tail  lights  flash 
synchronously  at  49.16  Hz,  beginning  at  about  one  second  after  gun  firing. 

This  results  in  loss  of  some  data,  but  the  delay  is  necessary  to  allow  the 
vehicle  to  clear  background  street  lighting.  The  flashes  are  recorded  by  three 
ballistic  cameras,  situated  as  shorn  in  figure  2.  The  position  data  is  obtained 
from  the  camera  records  by  the  method  outlined  in  reference  4. 

The  quality  of  the  trajectory  data  is  somewhat  lower  than  expected  from 
experience  with  the  previous  double  and  single  flashing  light  trials.  The 
estimated  r.m.s.  noise  levels  in  the  various  parameters  are  given  in  Table  2. 

The  main  cause  of  the  trouble  may  be  an  erratic  flash  rate.  Some  of  the  tail 
and  nose  flashes  were  inexplicably  delayed  and  flashes  were  not  simultaneous  for 
about  10  points  out  of  a total  of  170,  and  this  casts  doubt  on  the  general 
reliability  of  the  flash  triggering  unit. 

The  flight  history  of  the  vehicle  is  given  in  figures  3,  4 and  5.  The  record 
available  is  that  part  of  the  trajectory  for  which  a three  caaMra  solution 
is  possible.  In  any  case  the  incidence  oscillation  has  damped  to  a very  low 
amplitude  by  the  time  the  vehicle  moves  out  of  view  of  camera  one,  and  so  the 
data  following  this  are  of  no  interest.  The  position  and  velocity  histories 
are  shown  in  figure  3.  The  velocity  components  are  derived  from  the  position 
data  by  siiq>le  numerical  differentiation  techniques  outlined  in  reference  5. 

The  original  position  and  velocity  data  have  r.m.s.  noise  levels  which  are 
negligible  relative  to  the  scales  used  in  figure  3.  Hence  no  attempt  has  been 
made  to  show  the  individual  data  points.  The  differences  between  the  velocities 
obtained  as  model  output  from  the  parameter  estimation  and  the  measured  velocities 
are  also  negligible  relative  to  the  scales  of  figure  3.  Therefore  the  curves  in 
figure  3 represent  both  measured  data  and  model  outputs. 

Figures  4 and  5 show  the  angular  data  from  the  experiment.  Figure  4 gives 
a record  of  the  elevation  and  azimuth  of  the  longitudinal  axis  of  the  vehicle  and 
the  incidence  components.  Figure  5 shows  the  behaviour  of  the  vehicle  in  the 
incidence  plane.  The  measured  data  are  shown  as  points  on  these  two  figures  and 
the  scatter  is  appreciable.  The  corresponding  outputs  from  the  mathematical 
model  are  represented  on  the  fig^ires  by  the  solid  curves.  These  two  figures 
provide  an  indication  of  the  goodness  of  the  fit  obtained  with  the  parameter 
estimation  algorithm. 

An  outline  of  the  method  used  to  obtain  the  data  presented  in  these  figures 
is  included  here  for  completeness.  The  results  are  derived  from  simple 
geometrical  considerations.  The  basic  data  points  (x.j.,  y^,  z.j.)  and  (Xj^j,  yj^,  Zj^) 

are  obtained  by  triangulation  from  the  three  camera  records  using  a least 
squares  solution.  The  trajectory  (x,  y,  z)  of  the  vehicle  centre  of  gravity  is 
derived  from  these  data  using  the  relation 
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where  is  the  distance  of  the  vehicle  centre  of  gravity  aft  of  the  flashing 

light  at  the  nose,  expressed  as  a percentage  of  the  total  distance  between  the 
lights,  Ax  » - x^,  and  so  on.  The  velocity  (x,  y,  i),  relative  to  the  air, 

is  derived  from  the  position  data  using  the  methods  outlined  in  Chapter  V, 

Section  8 of  reference  5 for  differentiating  empirical  data  and  allowing  approx- 
imately for  wind.  The  attitude  angles  of  yaw  and  pitch  (or  azimuth  and 
elevation)  respectively  are  defined  by  the  relations 

tan  ^ = Ay /Ax 

tan  e = -Az/l(Ax)*  + (Ay)^!**  (2) 

and  the  distance  between  lights  is  given  by  the  relation, 

= (Ax)^  + (Ay)*  (Az)* . 

The  distance  ^ is  a good  measure  of  the  reliability  of  the  data  points  and  of  the 
general  r.m.s.  noise  level  in  the  position  data.  The  values  of  k,  Jr,  i,  ^ and  0 
are  used  as  data  for  the  parameter  estimation  technique.  The  values  of  incidence 
in  non-rolling  body  axes  can  be  derived  from  the  relations 


u 

cos  0 cos  yj)  cos  9 sin  \li  -sin  6 

i 

V 

-sin  ^ cos  ^ 0 

w 

sin  6 cos  (fr  sin  6 sin  i/r  cos  B 

i 

tan  f = w/u 

V 

tan  = v/u 

where  u,  v and  w are  the  components  of  vehicle  velocity  in  rotating  body  axes. 


3.  PARAMETER  ESTIMATION 

This  section  provides  a definition  of  the  mathematical  model  used  in  the  data 
analysis  and  outlines  the  method  of  data  analysis. 

3.1  Mathematical  model 

The  equations  of  motion  for  a rigid  body  with  six  degrees  of  freedom  in 
non-rolling  body  axes  are  given  in  reference  6 in  the  form, 

u ♦ qw  - rv  “ X/m  - g sin  0, 

■fr  ru  - Pw  * Y/m, 

w ♦ Pv  - qu  = Z/m  + g cos  0, 

® ■ q.  4 

^ ■ r/cos  6, 

4 - Pr(l  - I^/I)  • M/I, 

t * Pq(l  “ N/I,  (4) 
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where 


V =•  velocity  of  vehicle  in  rotating  body  axes. 


angular  velocity  of  vehicle  in  body  axes. 


angles  of  azimuth  and  elevation  respectively  of  vehicle  axis. 


« aerodynamic  forces , 


M » aerodynamic  moments, 

, N . 

P = -r  tan  B « residual  roll  rate  of  axes, 
m ’o  vehicle  mass,  and 

I^,  I » moments  of  inertia  of  vehicle  in  roll  and  pitch,  respectively. 

The  non-rolling  body  axes  have  the  OX  axis  forward  along  the  longitudinal 
body  axis,  the  OY  axis  in  a horizontal  plane  to  the  right,  and  the  OZ  axis 
completes  the  right  handed  set.  In  order  to  maintain  the  OY  axis  in  a 
horizontal  plane  a small  residual  roll  rate  must  be  maintained  for  the  axes 
system,  this  axes  roll  rate  is  called  P,  in  contrast  to  the  body  roll  rate, 
p,  which  does  not  appear  in  the  equations  because  of  the  assumption  that  the 
roll  equation  is  uncoupled  from  the  other  equations.  Since  there  is  no  roll 
information  available,  the  body  roll  rat*  and  roll  position  equations  are 
omitted  and  it  is  assumed  that  the  roll  dependence  of  the  aerodynamic  forces 
and  moments  is  negligible.  Ninety  per  cent  of  the  data  is  taken  for 
incidence  values  less  than  six  degrees  in  the  particular  example  treated 
here  so  that  this  assumption  is  unlikely  to  introduce  any  errors  larger  than 
the  r.m.s.  values.  Similarly,  althou^  the  first  non-linear  terms  in 
incidence  variation  of  restoring  moment  and  normal  force  are  included,  we 
cannot  expect  to  determine  them  reliably  because  there  is  insufficient  high 
incidence  data. 

The  expressions  for  the  aerodynamic  forces  and  moments  are, 

X • QSC„ 


Y * QS(C,  ♦ C_  tan*  f)  v/u 
Z - QS(C,  ♦ C,  tan*  f)  w/u 

Zf 


I 
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where 


M 

N 


QSd(  (C 
QSdl  (C_. 

"f 


tan*  f)  w/u  ♦ C (qd/2V)l 

B 

q 

tan*  f)  v/u  + C (rd/2V)l 

"q 


(5) 


tan*  f 


C tan  f C tan*  f = 

C_  tan  f ♦ C_  tan*  f = 

Zf  Zf* 


= (v*  + w*)/u* 

aerodynamic  pitching  moment  coefficient, 

aerodynamic  normal  force  coefficient, 

aerodynamic  axial  force  coefficient, 
aerodynamic  pitch  damping  derivative. 


Q = dynamic  pressure, 

S = force  and  moment  reference  area  and 
d = body  diameter. 


It  should  be  noted  that  the  Mach  number  is  less  than  0.4  so  that  the  flow  is 
incompressible  and  there  is  no  variation  of  any  of  the  coefficients  C^,  etc. 

with  Mach  number.  The  terms  C . C , C-  , C-  are  derived  from  the 

f ”t*  h 

Taylor  series  expansions  for  the  static  aerodynamic  pitching  moment  and 
normal  force.  The  vehicle  flies  at  low  subsonic  Mach  niimbers,  throughout 
the  flight,  and  we  have  assumed  that  there  is  no  roll  dependence  of  aerodyn- 
amic terms.  Therefore  the  pitching  moment  and  normal  force  coefficients 
are  functions  of  tan  f only.  As  the  majority  of  the  data  points  show  only 
low  values  of  incidence  (tan  f < 0.1)  we  can  use  Taylor  series  expansions 
for  normal  force  and  pitching  moment.  Allowing  for  the  fact  that  both  are 
odd  functions  of  f , we  obtain  the  expansions  used  above. 


C 


Z 


C_  tan  f ♦ C_  tan*  f ♦ . . . 


for  pitching  moment  and  normal  force  coefficients,  respectively 

Then  ^,0  are  obtained  directly,  as  model  outputs  and  the  velocities 
k,  f,  t are  obtained  from  the  simple  axis  transformation 


k 

cos  0 cos  ^ 

-sin  ^ 

sin  0 cos 

9 

m 

cos  0 sin  ^ 

cos 

sin  0 sin  4> 

t 

-sin  0 

0 

cos  0 

(6) 
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which  is  the  inverse  of  equation  (3).  These  five  quantities  can  then  be 
compared  with  values  for  the  same  variables  which  have  been  obtained  from 
the  trajectory  measurements.  The  initial  values  f t , q , r^ 

which  are  required  to  start  the  numerical  integration  of  equations  (4) 

are  parameters  which  must  be  determined  by  the  parameter  estimation  process, 

in  addition  to  values  for  the  aerodynamic  parameters  C-  , C-  , C , C 

X Zj3  0^3 

and  C The  parameters  are  referred  to  as  pi  , pj , ...  pi  3 in  the  rest  of 

”q 

this  report  and  are  taken  in  the  order  given  above,  where  pi  = x^,  and  so  on 
up  to  Pi  3 = C . Table  3 can  be  used  as  a guide  to  the  parameters 

"q 

represented  by  each  p^.  The  above  equations  and  parameters  thus  define  the 
whole  mathematical  model  used  in  the  parameter  estimation  process. 

3.2  Normal  equations 

The  technique  for  adjusting  initial  guesses  for  parameter  values  is 
essentially  a modified  Newton-Raphson  technique  as  discussed  in  reference  7, 
but  it  is  outlined  here  for  completeness.  If  there  are  measured  values 


■ij  = ®j(Pi  . Pi . 


Pn*  ^i^  ^ 


1.  2 ...  I,  j 


1 , ...  J 


of  J flight  variables  at  I times  t^,  where  IJ  > N,  the  parameters  pi  , pj , ... 
P|^  can  be  obtained  by  a recursive  solution  of  a set  of  normal  equations  which 
are  derived  below. 


Assuming  that  we  have  a current  set  of  estimates  p 


(s) 


1.  2 


N 


for  the  parameter  values,  then  by  using  the  mathematical  model  we  can  estimate 
a set  of  measurements  by  numerical  integration  of  the  model  equations,  and 
they  are. 


m.  . 
ij 


(s) 


Oj  (P«  . Pi  , . • • . tj)  • 

(s) 


This  set  can  be  used  to  form  a vector  of  residuals  R''  ''  where  each  element 


(s) 


where 


•ij  k ■ 1-  2.  •••  “ ■ k 


k « (i  - 1)  J ♦ j 


(7) 


Then  the  output  error  method  of  parameter  estimation  requires  the  the 
quantity 


i(s) 


l(s)T  ^ ^(s) 


(8) 


I 


is  minimised,  where  N is  a wei^ting  matrix.  The  matrix  N is  a K x K 
diagonal  matrix,  and  the  values  on  the  diagonal  are  proportional  to  the 
accuracy  of  each  measurement. 


WRE-TN-1719  (WRSD) 


- 7 - 


Suppose  that  new  parameter  values  can  be  represented  by  the  relations 


- (s+1)  ^ D ♦ 6d 


n' 


(9) 


so  that  the  new  set  of  estimates  for  the  measurements  can  be  written  in  the 
form 


m, 


(s+1)  _ 


”k 


(s) 


n 

E 

n=l 


ini 


(s) 


6p^  + 0l  (SPn)M  , 


(10) 


using  the  Taylor  series  expansion,  where  k is  defined  in  equation  (7), 
Now  the  sum  of  squares  of  residuals  is  given  by 

u(s+l)  , ,^(s+l)T  ^ ^(s+1) 


(11) 


where  " ™k  ” method  is  to  minimise 

The  minimum  occurs  when  (5pj^)  = 0 for  all  n.  The  Taylor  series 

expansion  from  equation  (10)  can  be  written  in  the  form 


M 


(s+1) 


(12) 


where  is  a matrix  of  partial  derivatives,  such  that  dj^^  = 

is  a vector  of  model  outputs  and  E is  a matrix  of  unknown  parameter 

adjustments  8p^,  Since  = M - where  M is  the  vector  of  measured 

values  the  above  equation  for  the  Taylor  series  expansion  can  be  rewritten 
in  the  form 


r(s.1)  . r(s)  . p(s),^ 


(13) 


,(s+l) 


Then  using  this  to  substitute  in  the  equation  (12)  for  U''  ■* , we  find  that 

u(s^l)  . W(R^®^  - D^*^E) 

and  hence  it  follows  from  the  conditions  for  a minimum, 

3y(»*l)/3(Jp^)  . 


0 for  all  n 
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-2  W + 2 W D^®^  E = 0. 


(14) 


Therefore,  if  we  represent  the  positive  definite  matrix  D^®^^  K D^®^ , by  "if, 
the  normal  equations  have  the  form 


'I'  E * W R^®^  , 


(15) 


which  gives  the  solution 


E = Sk  * D^®^^  W R^®^ 


(16) 


thus  yielding  values  for  the  parameter  increments  6p^. 

Some  statistics  are  available  on  the  accuracy  and  reliability  of  the  up- 
dated estimates  of  the  parameter  values.  The  following  results  are  quoted 
from  reference  7.  An  estimate  of  the  r.m.s.  value  a,  of  the  residuals,  or 

the  value  of  the  sum  of  squares  of  the  residuals  using  the  updated 

parameter  estimates  is  given  by 


o*  = o*  - E^  W R^®^ 

o " 


(17) 


where  o^  is  r.m.s.  of  the  residuals  of  the  observations  using  uncorrected 
parameters.  The  covariance  matrix  of  the  parameters  is  given  by 


cov  (E)  = o*  4'"'  , 


and  this  is  particularly  useful  in  determining  when  parameter  values  are 
correlated,  because  in  such  cases  estimated  values  of  the  correlated 
parameters  are  unreliable.  Finally  the  accuracy  of  each  parameter  estimate 
can  be  deteimined  from  the  standard  deviation  or  r.m.s.  error  level  <’(Pjj) 

which  is  given  by 


-*(P„) 


nn 


where  is  the  n^^  diagonal  element  of  . 
nn 


(18) 


3.3  Partial  derivatives 

An  important  requirement  of  the  parameter  estimation  technique  outlined 
above  is  to  estimate  the  partial  derivatives  of  the  estimates  of  the  observed 
values,  with  respect  to  the  parameters.  There  are  alternative  methods  of 
approach,  but  the  most  effective  presently  available  seems  to  be  to  use 
partial  differentiation  of  equations  (4)  and  (S)  which  define  the  mathematical 
model  in  Section  3.1.  This  approach  yields  the  following  simultaneous 
differential  equations,  for  j ■ 1,  2,  ...  13, 


9 
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d . 3w  ^ dv  Qs  * a 


J ‘ J 


Jt  (Is)  =-■■5?.  -“5^.  •'■55.  *“5^.  ‘^I(P.  -p..  tan’OfJ, 

Fj  Fj  Fj  Fj  Fj  Fj 

* (*9j  ^ ®10j  ''  2p.o  ^2  Cv  + W |^_)1 


^ ^Qs 


iy 


(19) 


d ,3w  . „ 3v  dp  du  ^ QS  f , ...  j F^ 

dt  d^.  - di.  ■"  ‘I  di.  5?.  ^ I fP’  * P">  di. 


1 ‘ 1 


* ^*9j  ^ *10j  ^ 2P'®  7 ip . ^ “ ip.^^  - g sin  0 


d0 


dt  ^dp.^  " 3B. 

J 


dr 


d0 


f-t  (^  ) = sec  0 ^ + r sec  0 tan  0 g-^ 

j J 


Pj 


Jt  fjJ.)  = %.  " ^ {(p*-  " P‘^  ii. 

*^1  J J J 

^ ^®llj  * *12j  * 2p,2  Jlr  (V  + w g^_) 


+ *4  Pi3  d + 5s  qd 


d dr  . 
dt 


= -d  - I/I)(P  q Tp^)  > 1-(P..  ^ P.4  tan^  f) 


- (*llj  ^ *12j  n V - 2p,2  Js  (V  + w |j_) 


♦ *S  Pi3  d Ij  + »S  dr 


I 

I 
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where  dP/dp^  « - tan  9 - r sec*  036 /3p^.  The  partial  derivatives 

for  ^ and  0 cone  directly  fron  this  set  of  differential  equations.  The 
partial  derivatives  for  x,  y and  i must  be  derived  fron  the  axes  transfom- 
ation,  given  in  Section  3.1, 


X 

cos  0 cos  ^ 

-sin  0 

sin  0 

cos  ^ 

u 

9 

cos  6 sin  ^ 

cos  ^ 

sin  0 

sin  ^ 

v 

i 

-sin  0 

0 

cos  6 

w 

Partial  differentiation  of  this  relation  yields 


di. 

dpj 


1^  cos  6 cos  sin  ® ^ ^ 


♦ i cos  w 5“ 
aPj 


cos  0 sin  ^ 5^  cos  ^ sin  0 sin  ^ * k ^ 


X • (3® 

+ i sin  w 5“ 

3p^ 


(20) 


3i 

3p 


3 


- sin  ® cos  0 - (u  cos  0 + w sin  0)  ^ 


3P^ 


The  set  of  simultaneous,  ordinary  differential  equations  in  partial 
derivatives  given  in  equation  (19)  with  the  equations  of  motion  quoted  in 
equation  (4)  comprise  the  complete  mathematical  model  and  form  the  complete 
set  of  98  first  order  simultaneous  differential  equations  (for  13  parameters) 
which  must  be  solved  on  each  iteration  of  the  parameter  estimation  method 
so  that  98  initial  values  are  required  to  start  the  solution.  The  first 
six  parameters  represent  the  initial  values  0^,  q^,  r^  in 

that  order.  The  inverse  transform  to  that  stated  above  is  required  to 
find  (u^,  v^,  w^)  and  the  inverse  is  obtained  simply  by  transposing  the 

matrix  so  that 


u 


cos  0 cos 

cos 

0 

sin 

0 

-sin 

o 

cos 

0 

'^o 

0 

sin  0 cos 
o 

sin 

0 

o 

sin 

0 

-sin  Oj 
0 


'<• 

o 

i 

o 

(21) 


The  initial  values  for  some  of  the  partial  derivatives  are  rather  more 
complex.  First  the  simple  ones  are 


'Jitfttt  ■ 
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where  6^^^  is  the  Kronecker  delta  symbol;  1 when  i = j and  0 otherwise. 

The  other  initial  conditions  must  be  obtained  through  the  axes  transformation 
and  are  as  follows. 


S’?, I ^^0  • ®o  ^o* 


I « cos  6^  sin  I ■ -sin  0^, 

3p2l  t=o  ° ^ t=0  ® 


51  ■ V cos  0 , 

^P4  •_«  O O 


Sps' 


^Pjl  t.( 


* 0 for  j > 6, 


t»0 


■sin 


* cos  tl/  , 
^p»'  t=0  ° 


t.o  ■ "■ 


^ I * COS  sin  I =0, 

*P*'  t»0  ° ° ° ° ®P*  t-0 


dp . 

*'y  t 


0 for  j > 6, 


hi 

^ n - sin  6^  cos  ^ I = sin  0 sin  ^ , I?  1 = cos  0, 

t=0  o ’^o’  op2't_Q  o o’  opal  { 


51  1 = V sin  0 , 

3p4'  o o’ 


^5 1 " “o 


= u and 


5^  I 

t=0 


0 for  j > 6. 


3.4  Algorithm  implementation 

The  algorithm  proceeds  as  follows.  We  wish  to  set  up  and  solve  the 
normal  equations  (IS),  which  are  a set  of  N linear  simultaneous  equations  for 
the  unknown  perturbations  in  the  parameter  values  p^,  j = 1,  2,  ...  N.  The 

normal  equations  (15)  can  be  rewritten  in  the  form 

A 3c  ■ b 

where  x is  the  column  vector  of  unknown  parameter  perturbations,  A is  the 
matrix  of  partial  derivatives  and  ^ is  a column  vector  of  residuals.  The 
matrix.  A,  has  members 


^ 5 ^ 5 a_  (s)  3 (s) 

i-l  j-1  k-1  *-l 


and  the  vector,  b,  has  members 


I 5 3 (s) 

EL*^  '•!)  - "j 

i-i  j-1 
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where  n,  m « 1,  2 


The  measurenents 


N,  where  N is  the  number  of  parameters  being  sought, 
represent  the  measured  values  of  i,  i,  and  6 at  time 

p,  for  i » 1,  2,  ...  I and  so  on  up  to  m^^  = 0(tj), 

The  weights  w^  for  j * 1,  2,  ...  5 are  proportional 

to  the  accuracies  in  the  measured  variables  i.,  t,  ^ and  0.  The  values 

(s) 


t. 

I 


so  that  m^j 


for  i * 1,  2, 


'/dp  can  be  obtained  from  the  numerical 
n 


for  the  partial  derivatives  dm^^^ 

solutions  of  the  mathematical  model  equations  (19),  for  the  partial 
derivatives  and  by  using  the  transformation  equation  (20).  The  model  outputs 


(s) 


used  to  evaluate  the  vector  of  residuals,  b,  can  be  obtained  from  the 


m.  . 
ij 

numerical  solutions  of  the  mathematical  model  equations  (4)  and  the  axes 
transformations  equation  (6) . Then  the  normal  equations  can  be  solved  for 
the  perturbations  to  the  parameter  values  and  the  next  iteration  of  the 
algorithm  begun. 

The  normal  equations  discussed  above  simply  tell  us  how  to  inq}rove  our 
estimates  of  the  parameter  values.  Two  major  problems  remain.  First,  how 
to  choose  first  values  for  the  parameters,  and  second,  how  to  control  the 
algorithm  so  as  to  encourage  convergence  rather  than  divergence.  The 
problem  of  choosing  initial  values  is  not  very  difficult.  The  initial  values 
of  velocity,  angular  displacements  and  angular  velocities  can  be  estimated 
reasonably  well  from  the  data.  Apart  from  this,  only  C is  really  required; 

"f 

the  remaining  values  can  be  started  at  zero.  The  pitching  moment  derivative 
C can  be  easily  estimated  from  the  relation 


w ’ = QSd  C /I 

"f 

where  cj  is  the  frequency  in  radians  per  second  of  the  pitching  oscillations 
shown  in  figures  4 and  5.  Methods  for  obtaining  convergence  are  difficult 

to  express  in  simple  general  rules,  since  they  will  vary  with  the  problem 
being  treated.  However  a basic  technique  is  to  avoid  trying  to  find  values 
for  all  parameters  at  once.  A gener  1 approach  is  to  arrange  the  parameters 
in  order  of  importance  of  their  effeci.  on  the  motion  or  model  outputs. 

Then  only  the  first  two  or  three  should  be  allowed  to  vary  and  convergence 
obtained  for  that  situation.  When  convergence  has  been  reached  two  or  three 
of  the  other  parameters  should  also  be  allowed  to  vary  until  convergence  is 
achieved  once  again.  The  process  should  be  allowed  to  proceed  in  this  way 
until  all  parameters  have  been  determined  and  the  solution  optimised. 


4.  RESULTS 

The  paraawter  values  obtained  frra  the  pararoter  estimation  analysis  of  the 
trials  results  are  given  in  Table  3,  together  with  the  estimated  root  mean 
square  errors  in  each.  These  values  were  obtained  using  150  points  out  of  the 
total  of  157  which  were  finally  available,  beginning  at  the  seventh  flash. 

A weight  of  0.1  was  used  in  matching  against  the  velocity  components  relative  to 
a weight  of  unity  for  the  angular  data.  This  can  be  related  to  the  r.m.s. 
noise  levels  given  in  Table  2 which  show  in  broad  terms  that  noise  levels  are  an 
order  greater,  nimwrlcally,  in  velocity  components,  than  in  attitude  angles, 
althou^  much  less  in  percentage  terms. 

The  accuracy  of  the  fit  obtained  in  terms  of  measured  data  and  model  outputs 
is  apparent  from  figures  3,  4 and  5.  Figure  3 shows  velocities  and  on  the 
scales  used  the  differences  between  measured  values  and  model  outputs  are 
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negligible.  Figures  4 and  5 show  angular  data  both  attitude  angles  and 
incidence  angles.  Here  the  noise  in  the  data  is  appreciable  and  the  scatter  of 
the  data  points  compared  with  the  amplitude  of  the  oscillation  provides  a good 
idea  of  how  well  the  parameter  estimation  method  has  matched  the  data.  The 
figures  also  show  that  substantial  improvements  can  be  expected  in  the  accuracy 
of  derived  aerodynamic  coefficients,  if  the  noise  levels  in  the  angular  data 
are  nearer  the  expected  value  of  one  degree,  rather  than  three  or  four  degrees. 

Table  3 shows  that  initial  values  of  state  variables  are  well  determined. 
However  the  quantities  of  main  interest  are  the  aerodynamic  coeffi'-.ients. 

The  drag  coefficient  is  well  determined  with  an  r.m.s.  error  of  only  2.2%,  as 
is  also  the  pitching  moment  derivative  at  zero  incidence,  C , with  an  r.m.s. 

mj. 

error  of  0.9%.  The  normal  force  derivative  near  zero  incidence  and  the  pitch 
damping  derivative  are  less  well  determined,  with  r.m.s.  errors  around  7%. 

The  derivation  of  normal  force  depends  mainly  on  velocity  measurements  and  the 
insensitivity  of  the  changes  in  velocity  to  this  parameter  is  reflected  in  the 
low  accuracy  of  the  normal  force  derivative.  Since  90%  of  the  data  was  taken 
at  incidence  amplitudes  less  than  six  degrees  the  information  on  non-linear 
normal  force  and  pitching  moment  terms,  C.,  and  C , was  minimal.  In  fact, 

mj.3’ 

attempts  to  determine  the  non-linear  terms  produced  values  with  r.m.s.  errors 
of  between  20  and  50%  while  simultaneously  doubling  the  uncertainty  in  the 
linear  terms.  Therefore  the  attempts  were  abandoned. 

Some  effort  was  made  to  verify  the  results  by  comparison  with  wind  tunnel 
measurements  given  in  reference  3.  The  principal  wind  tunnel  results  are 
plotted  in  figures  6 and  7 which  show  pitching  moment  coefficient  and  axial  and 
normal  force  coefficients  respectively.  Unfortunately  there  were  some 
differences  between  the  flight  model,  shown  in  figure  1 and  the  model  used  in 
wind  tunnel  measurements.  The  first  and  most  obvious  is  the  cylindrical  cap 
which  houses  the  rear  flashing  light.  However,  it  is  estimated  that  this  will 
not  produce  appreciable  effects  on  the  aerodynamic  coefficients.  The  second 
difference  which  is  more  important  is  that  the  fins  on  the  flight  model  are  only 
0.018  calibres  thick,  whereas  the  fins  on  the  wind  tunnel  model  are  0.031 
calibres  thick.  Consequently  the  chamfer  on  the  leading  edge  of  the  flight  model 
reaches  only  about  halfway  to  the  junction  of  the  boattail  and  cylindrical  tail 
section  of  the  body,  while  the  chamfer  on  the  wind  tunnel  model  reaches  the 
junction  exactly.  This  raises  the  possibility  of  separation  from  the  leading 
edge  of  the  fins  due  to  the  unfavourable  pressure  gradients  produced,  and  a 
consequent  loss  in  efficiency  of  the  tail  of  the  free  flight  model.  Two  other 
differences  between  wind  tunnel  tests  and  the  flight  test,  which  would  probably 
affect  drag  only  are  the  Reynolds  number,  and  the  roughness  band  which  was  used 
on  the  wind  tunnel  model  to  induce  a turbulent  boundary  layer.  The  Reynolds 
number  based  on  body  diameter  was  2.94  x 10^  for  the  wind  tunnel  tests  and  varied 
from  1.437  x 10*  down  to  1.317  x 10*  in  the  flight  tests. 

Let  us  look  first  at  the  axial  force  coefficient  C^.  The  value  obtained 

from  wind  tunnel  measurements  is  " 0.123  and  the  value  obtained  from  this 

free  flight  experiment  is  (C^^jpp  » 0.100,  a value  which  is  further  siqpported  by 

the  measurement  of  0.097  obtained  from  the  previous  trial  reported  in  reference  1. 
As  far  as  axial  force  is  concerned  the  main  difference  between  the  wind  tunnel 
measurements  and  the  free  flight  trials  is  the  Reynolds  number.  The  skin 
friction  contribution  to  the  drag  on  a flat  plate  varies  with  Reynolds  number, 
according  to  the  relation 

Cj«(log,o  Re)"^-^^ 
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Since  the  skin  friction  contributes  about  8S%  of  the  axial  force  on  the  model, 
the  corrected  value  of  the  free  flight  measurement  for  a Reynolds  number 
equivalent  to  the  wind  tunnel  value  is  ■ 0.126.  Allowances  for  the 

effects  of  thinner  fins  and  absence  of  a roughness  band  on  the  free  flight  model 

would  increase  this  estimate  slightly.  However,  the  differences  between  the 

value  corrected  for  Reynolds  number  and  the  value  measured  in  the  wind  tunnel 

would  remain  within  the  experimental  error  limits. 

The  normal  force  derivative  measured  in  free  flight  was  (C_  » -3.5, 

Ft* 

while  wind  tunnel  measurements  gave  (C^  - -4.2  a difference  of  0.7. 

However,  the  r.m.s.  error  in  the  flight  measurement  is  0.24,  so  that  the 
difference  is  less  than  three  standard  deviations,  assuming  that  the  error 
distribution  is  normal.  While  such  a discrepancy  is  by  no  means  probable  it 
is  not,  statistically,  highly  significant.  On  the  other  hand  flight  measurements 
gave  (C^)pp  “ -3.8.  The  discrepancy  here  is  0.5  compared  with  an  r.m.s.  error 

of  0.03,  so  that  the  difference  is  clearly  significant.  However,  since  the 
moment  arm  for  the  tail  lift  is  about  twice  that  for  the  nose,  relatively  small 
variations  in  tail  efficiency  result  in  relatively  large  changes  in  static 
stability,  and  so  it  seems  likely  that  this  discrepancy  can  be  accounted  for  by 
some  small  additional  separation  on  the  fins  of  the  free  flight  model  due  to 
the  slight  differences  in  the  leading  edge  chamfer.  In  fact  only  seven  per  cent 
loss  of  efficiency  of  the  tail  will  lead  to  the  values  of  C obtained  from  the 

“r 

free  flight  experiment.  This  results  in  a loss  of  only  3.5%  of  total  lift. 
Examination  of  figure  1 shows  that  lift  could  easily  be  lost  from  seven  per  cent 
of  the  fin  area  due  to  separations  arising  from  the  fact  that  the  leading  edge 
chamfer  falls  short  of  the  end  of  the  boat  tail.  Such  a result  is  of  course 
also  consistent  with  the  value  obtained  for  C_  , although  the  r.m.s.  errors  in 

C-  are  so  large  that  nothing  definite  can  be  said  about  values  obtained  for 


normal  force. 

Thus,  as  far  as  can  be  ascertained,  the  free  flight  experimental  results  are 
consistent  with  wind  tunnel  measurements.  Although  there  are  some  discrepancies, 
they  can  be  accounted  for  by  differences  between  the  two  models.  In  addition, 
the  axial  force  coefficient  measurement  compares  well  with  the  value  of  0.095 
derived  from  a previous  flight (ref. 1)  during  which  the  model  achieved  no 
significant  incidence. 


5.  CONCLUSION 

A simple,  cheap  and  convenient  technique  for  flight  testing  of  missile  shapes 
has  been  presented.  The  advantages  of  the  technique  arise  from  the  large  amount 
of  information  which  can  be  derived  from  relatively  few  measurements.  The 
only  on  board  facilities  required  are  two  flashing  lights  and  a small  flash 
initiator  unit.  The  vehicle  is  launched  from  a compressed  air  gun,  and  the 
light  flashes  are  recorded  by  three  ballistic  cameras.  Values  for  all  the 
important  aerodynamic  derivatives  can  be  derived  from  the  experisttntal  data, 
using  a parameter  estimation  technique.  The  steps  in  gathering  and  analysing 
the  data  are  all  si]q)le  and  straightforward  although  -lOme,  such  as  reading  the 
camera  plates  are  soswwhat  onerous  and  exacting. 

Comparing  the  free  flight  results  with  wind  tunnel  measurements  shows  some 
discrepancies.  However,  the  difference  between  wind  tunnel  values  and  free 
flight  values  for  the  aerodynamic  derivatives  can  be  explained  by  differences 
between  the  two  models  used,  and  by  Reynolds  number  effects. 
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At  this  stage  of  development  of  the  flashing  light  technique  it  is  possible 
to  measure  all  the  significant  aerodynamic  derivative  values  near  zero  incidence, 
that  is,  all  the  significant  linear  aerodynamic  terms.  The  next  step  is  to 
take  some  account  of  the  non-linear  aerodynamic  effects  which  become  significant 
at  higher  values  of  incidence.  Non-linear  aerodynamic  effects  are  always 
dependent  on  the  orientation  of  the  vehicle  relative  to  the  incidence  plane. 
Hence,  in  order  to  extend  the  usefulness  of  the  technique,  a means  of  measuring 
the  roll  orientation  of  the  vehicle  needs  to  be  developed.  This  may  require 
the  inclusion  of  telemetry  in  the  system,  although  this  should  be  avoided  if  at 
all  possible  because  both  the  free  flight  vehicle  and  the  data  collection  and 
analysis  will  then  become  much  more  complicated  and  expensive.  However,  if  the 
method  is  to  be  generally  useful  some  means  of  measuring  roll  orientation  of  the 
vehicle  will  have  to  be  found. 
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NOTATION 

skin  friction  coefficient 
axial  force  coefficient 
linear  noraal  force  coefficient 

non  linear  normal  force  coefficient 

linear  restoring  moment  coefficient 

non  linear  restoring  moment  coefficient 

pitch  dangling  coefficient 

body  diameter 
gravitational  acceleration 
moment  of  inertia  in  roll 
moment  of  inertia  in  pitch 

total  aerodynamic  moments 

vehicle  length  - distance  between  lights 

distance  of  centre  of  gravity  aft  of  nose  light,  as  a 
fraction  of  total  distance  between  lights 

vehicle  mass 

measured  values 

residual  roll  rate  of  body  axes 

angular  velocity  vector  of  body 

Reynolds  nianber 
parameters  to  be  estimated 
changes  to  parameters 
dynamic  pressure 

residual  errors  in  estiswted  measurements 
reference  area  for  aerodynamic  forces  and  moments 
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a 

♦ 


OJ 

superscripts 


(s) 

subscripts 

o 

N 

T 


time 


vehicle  velocity  vector  in  rotating  body  axes 


total  aerodynamic  forces 


vehicle  position  in  range  axes 


differences  between  nose  and  tail  positions 

total  vehicle  incidence 
horizontal  incidence  component 
vertical  incidence  component 
elevation  angle  for  vehicle  attitude 
r.m.s.  error  levels 
matrix  in  normal  equations 
azimuth  angle  for  vehicle  attitude 
frequency  of  pitching  oscillation  (rad/s) 

differentiation  with  respect  to  time 
estimate  using  current  parameter  values 

initial  values 

values  pertaining  to  nose  light  trajectory 
values  pertaining  to  tail  light  trajectory 
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TABLE  1.  MISSILE  PHYSICAL  DATA 


Mass 

17.9  kg 

Roll  inertia 

0.046  kg  m* 

Pitch  inertia 

1.34  kg  m^ 

Diameter  (calibre) 

0.175  m 

Reference  Area 

0.02405  m^ 

Centre  of  gravity 

2.39  calibres  from  nose 

TABLE  2.  MEAN  R.M.S.  NOISE  LEVELS  ON  TRAJECTORY  DATA 


Variable 

Nose 

Tail 

CG 

X (m) 

0.013 

0.038 

- 

y (m) 

0.010 

0.048 

- 

z (m) 

0.018 

0.065 

- 

k (m/s) 

0.47 

1.23 

0.47 

9 (m/s) 

0.49 

1.11 

0.45 

i (m/s) 

0.32 

0.86 

0.40 

>1/  (rad) 

- 

- 

0.045 

e (rad) 

0.071 

Note:  The  length  of  the  vehicle  calculated  from  nose  and  tail 

trajectories  has  a mean,  1.077  m,  and  a standard  deviation, 
0.0280  m.  If  we  assume  that  nose  and  tail  trajectories 
contribute  equally  to  the  standard  deviation  then  the 
range  coordinates,  which  are  the  most  significant  in 
determining  vehicle  length,  must  have  an  r.m.s.  noise 
level  which  averages  0.0198  m.  Although  a little  large, 
this  is  reasonably  consistent  with  the  above  table. 


Figure  1.  Test  vehicle 
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Figure  2.  Range  layout 
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Figure  3 


(,.Siu)  Ail0013A  *''' 


(UJ)  3DNViSia  *-•< 


(UJ)  33NViSia  X 


Figure  3.  Trajectory  of  vehicle 


WRE-TN-1719  (WR6D) 
Figure  4 


Figure  4.  Attitude  and  incidence  histories 
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Figure  5 


Figure  5.  Incidence  behaviour 
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Figure  6 
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Figure  7 


Figure  7.  Axial  and  normal  force  coefficients 
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1 Summary 

'V This  report  describes  a new  method  for  obtaining  the  low  incidence, 
linear  aerodynamic  data  for  a missile  from  very  simple  free  flight  trials 
on  a gas  gun  range.  The  only  instrumentation  in  the  free  flight  model 
is  two  flashing  lights,  one  in  the  nose  and  one  in  the  tail.  The  images 
of  these  flashing  lights  are  recorded  by  three  ballistic  cameras  and  the 
positions  of  the  lights  can  then  be  calculated  very  accurately  by  a least 
squares  solution  of  the  resulting  triangulation  problem.  The  velocity 
of  the  centre  of  gravity  and  the  attitude  of  the  vehicle  can  be  estimated 
from  the  position  measurements  for  the  nose  and  the  tail.  A parameter 
estimation  technique  using  the  output  error  criterion  is  applied  to  the 
attitude  and  velocity  data  to  obtain  aerodynamic  terms.  The  method  has 
been  used  successfully  to  estimate  axial  force,  normal  force,  static 
stability  and  pitch  damping  moment  for  the  test  vehicle.*. 
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